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Abstract 

The probability density function of the positive occupation time of one-dimensional 
Brownian motion with two-valued drift is derived by using the Feynman-Kac formula and 
Laplace transforms. The result is applied to a simple, two-dimensional stochastic differential 
equation describing stochastically perturbed sliding motion of a discontinuous, piecewise- 
smooth vector field. Long time asymptotics of the density are also computed. 

1 Introduction 

The one- dimensional stochastic differential equation 

= I ^ ;J: > } dt+dw{t) , x(o) = xo , (i) 

where aL,(iR G K and W{t) is a standard Brownian motion, represents Brownian motion with 
two-valued drift. Equation ([T]) has arisen in a stochastic control problem [1], investigations into 
the effects of noise on sliding motion in piecewise-smooth systems [2], and general theoretical 
settings |3lll]. In the context of sliding motion it is assumed ai.aR > 0, such that solutions to 
(II]) rarely stray far from the origin, but in this paper such a restriction is not made. An explicit 
expression for the PDF (probability density function) of x{t) was first derived by Karatzas and 
Shreve in [5]. 

For any t > 0, let 

T = / X[o,oo)(a;(s)) c/s = meas {s G [0,t] I x(s) > 0} , (2) 
Jo 



*The author thanks Rachel Kuske for many helpful discussions regarding this work. 
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denote the positive occupation time of x{t). The purpose of this paper is to derive the PDF of r, 
call it p, i.e., 

Prob(r G f]) = / p{f;t;xo,aL,aR)df, (3) 
Jn 

for any measurable subset, Q C [0,t]. PDFs of occupation times of simple stochastic processes 
have found applications in control problems and mathematical finance [6]. 

The theorem below gives p for xq = 0. The result for xq may be written in terms of the 
result for a;o = by conditioning over the first passage time of x{t) to zero. Specifically, let 

/i(s,xo,a) = — ==e 2s , (4) 

denote the PDF for the first passage time, s, of dx = adt + dW{t), from a;(0) = xq to x(t) = 0. 
Then 



p{T;t;xo,aL 



^ 1 /(j a;o, aL)p(r; t - s; 0,0^,0^) ds + 5{t) /^°° Xq, Cl) ds 
' ^ I /o Xq, -aR)p{T - s;t-s;0, aL, qr) ds + 6{t - t) J^°° h{s, : 



Xq <0 

^0, ~(^r) ds , Xq > 

(5) 

where 6 is the Dirac-delta function and we set p(r; t; 0, ai, a,R) = whenever r ^ [0, t]. For both 
Xq < and Xq > 0, the first term of ([5]) corresponds to the case that x(t) first reaches zero at a 
time s <t, and the second term corresponds to x(t) not reaching zero by the time t. 

Theorem 1. 

For xo = and any t > 0, the PDF of the positive occupation time ^ of Brownian motion with 
two-valued drift is 



p{T;t;0,aL,aR) 



— -Q-^ «Re 2 erfc a^a 2 erfc y 



■n^rit - r) A/27r(t - r) yfhvr 



, \Pl{aL^aR) -{{aL^o.R)r-aj^tf I -{a^ + aR)Wr{t - t) 
H = e 2t erfc ' 



^TTt y V2t 

+ J^{t; t] ai, aR) + J^{t - r; t; aR, ai) , (6) 



where 



-a2 T -a2 z 

^" " 20F io v^(^ + r) +V2(. + r)l 1^ V^{^) J 

(7) 



In two special cases, (E]) reduces to known results. First, if a^, = = 0, then x{t) is 
simply regular Brownian motion. As first shown by Levy in this case the occupation time is 
distributed by the arc-sine distribution [HI [9] . Indeed, in this case only the first term of ([6]) is 
nonzero, and 

p(r;t; 0,0,0)= .} . (8) 
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Equation ([8]) provides a reasonable approximation to p{T;t]0,aL,aii) for small t, see Fig. [H 
because over short time frames the diffusion of a stochastic quantity dominates its drift. 

Secondly, if = — a/j, then x{t) is Brownian motion with constant drift. In this case 
simplifies to 

p(r; 0. -a. a) = (-i= e^*' - ^ erfc (^) ) (^^=j + ^ (^^^) ) 

(9) 

which was first derived by Akatori [lU]. The simplification is not straight-forward but may be 
demonstrated by using different expressions for Owen's T- function [TT] . 

The remainder of this paper is organized as follows. A proof of Theorem [1] is presented in 
[|2j The long time asymptotics of are described in ^ focusing on the case a^, a^j > 0. In §11 
the result is applied to a model of stochastically perturbed sliding motion. Finally ^ contains a 
brief discussion. 




Figure 1: The PDF ([6]) when ol = 2 and = 1. Note that the axes are scaled by t such that 
the domain is the unit interval at all times. For small t, p is well-approximated by the density of 
the arc-sine distribution, ([H]). For large t, p is approximately Gaussian since ai^a^ > 0, see §21 
When t = 10, on the given axes the Gaussian approximation has mean ^^^^^ = | and standard 
deviation ^ , ^ = -4= ^ 0.11, as shown. 

(aL+anjVt 3V10 ' 
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2 Proof of Theorem [T] 

For all t > 0, the stochastic differential equation ([T]) has a unique strong solution US] and 
therefore the expectation 

where A > 0, is well-defined. Equation ([6]) may be obtained from u because we have 

u{0,t,X)= [ e-^^p{T;t;0,aL,aR)dT . (11) 



Following a method applied to similar problems, see for instance O [TOl [T3], by the Feynman- 
Kac formula P |15], u is the unique, bounded, continuous solution to the initial value problem 



Integration of the PDF with respect to over an arbitrarily small neighbourhood of zero reveals 
that is also continuous. 

To solve f|T2|) . we take the Laplace transform 

POO 

U{xo,a,X)= / e-''^u{xo,t,X)dt , (13) 
Jo 

to produce the piecewise-linear ordinary differential equation 

_l ^ I Wxoxo + o-lUxo - all , xo < 

\ ^Uxoxo - anUxo - (a + X)U , xq > ' 

Via standard ODE methods, we can obtain an explicit expression for the unique bounded solution 
to f|T^ for which u and Ux^ are continuous at xq = 0, and from this arrive at 

a ( A/2(a + A) + a| - ) + (a + A) ( ^/2a + a\ - ) 

t/(0, a. A) = ^ ^ \ ^ . (15) 

a{a + A) ( y2{a + X) + a\ — + -^2a + a\ — a^] 

Our goal is to obtain p, which by (ITT!) and (fT3|l is related to [/(O, a. A) by 

U{0, a. A) = / / e-°*e-^Xr; t; 0, a^, ur) dr dt . (16) 
Jo Jo 

Equation 0161) is now used to derive p constructively. To make the task simply the evaluation of 
two inverse Laplace transforms, we first reverse the order of integration in (fT6l) . and let f = t — r 
and fi = a + X, to obtain 

POO POO 

t/(0,a,/i-a) = / / e-'^''e-''^p{T,v + T;0,aL,aR)dvdT . (17) 
Jo Jo 



Simplification is provided by noting that 

U (0, a,fi~a) = E{fi, a; ql, a^) + E:(a, /i; a^, a^) , (18) 

where 

^(/i, a; ax, ai?) = . ^ r- • (19) 

By (fT71) . with a httle care it follows that 

p(r; t; 0, a^, a^j) = V(r, t - r; a/,, a^j) + ^{t - r, r; a^j, a^) , (20) 

where 

/•oo /»oo 

a; ttL, cLji) = / / e~°''"e~^'^il){T,v]aL,aR) dv dr . (21) 
Jo io 



We now use fl2T|) to obtain ^/^ from f|T9|) . Using, in particular, £ ^ (^-^7=:^; j = * {k\/t) , 

the inverse Laplace transform of f|T9l) with respect to a is 

A/2/i + a| - / V2 

2/i I 0™ 



e ^''tp{T,v]aL,aR)dT = ^ e 2 { ^== - { \/2fi + a% - aR - ai 



X e^ 



^(^^-"™)^^erfc((v'^:74 ^) ^22) 

To perform the daunting inverse Laplace transform of fl22|) . we rewrite it as 
/"°° _ z4z / 1 

/ e ''''^{T,v;aL,aR)dT = e 2 [ ^=Gi{fi;aR) ~ G2{fJ,,v;aL,aR) 

Jo Vv^^ 



+ ^^^^^^^^^Gi{fi]aR)G2ifi,v;aL,aR) ) , (23) 



where 



G,(^;.h) = V^-o^ . (24) 
G2(fi, t>;ai,aH) = e(\/*+°I-'"«-'")'* erfc (|(^\/2,^ + aj, - - a^^ ^) , (25) 



SO that 



tp{T, v; ai, aR) = e 2" ( gi{T; aR) - g2{T, v; ai, aR) + ^^^^ 9i{'r; clr) * g2{T, v; ai, aR] 



(26) 

where gi and g2 denote the inverse Laplace transforms of Gi and G2 respectively, and * denotes 
convolution with respect to r. We have 



^i(r;a^) = ^-^erfc(^-^) , (27) 
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and by the identity £ ^ fQiVs+^)'^ erfc {^/s + ku = -^-^ J°° ze * (""^"^Je c/z (obtained by com- 
bining formulas in [IS]), we also have 

^2 r,t;;aL,aR = — — -j= + e 2(r+„) erfc , , • 28 

After an involved algebraic simplification omitted for brevity, the combination of fl20|) . fl26|) . fl27j) 
and fl25]) yields ([H]) with ([7]) as required. □ 



3 Long time asymptotics of p 



This section provides a statement of the long time asymptotics of ([6]), then a discussion of its 
derivation. 
For large t, 



p{T;t;0,aL,aR) ~ < 



^(r; ol, ctij) , < 0, a/j > ^ (29) 

Q(t - T;aR,aL) , ai > 0, < 

^ Q{t; aL, a^) + g{t - t; a^, a^) , ai <0, < 



where, 



gir; a„ a«) = e^ - a,{2a, + an) e^^^— )^ erfc ^zi^^^Ll^^ . (30) 

To obtain fl29|) . first note that it may be shown directly from ([6]), that for any fixed r > 0, 

p{T;t;0,aL,aR) ^ -I > ' 

as t — >■ 00. For brevity we omit a derivation of this statement which requires taking t — )■ 00 in 
(I7j) and simplifying the result via integral transformations. Moreover, for ql < 0, 

r „ < n 

^(r;a^,a«)dr= ' Z n ^ (32) 

which ensures (l29i) has unit area in the limit t — )■ 00. The term Q(t — t; or, a^) appears in ( 129|1 
by symmetry. 

When aL,aR > 0, for large t the first three terms of ([6]) may be neglected because their Li 
norms rapidly approach zero as t — t- 00. The fourth term of ([6]) approaches ^f2^i ^ 
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because here erfc(-) — )■ 2. To simplify ([7]), we consider the two terms in the integrand of ([7]) sep- 
arately. The first term is negligible for large t when aL-,ciR > 0. By letting 



be a scaled difference of r from tlie mean, ^^^l^^ , wliere (3 > —I since r G [0,t], and substituting 
u = we liave 

2 

r(r; t, a„ a,) ^ "^""^ t^"'^ , e^-^'^H erfc(.) A 



2v27r / aL^fl "fip 



(34) 

Equation fl34|) may be evaluated asymptotically [17]. When p < 0, this is achieved by expanding 
the integrand in a Taylor series about u = 0; when p > 0, we expand about the upper limit of 
the integral. For instance, when p < (the other case is similar and produces the same result), 
after integration we have 



T{T;t;aL,aR) , ^ e v t^+i ; + o I - 1 



(35) 



Only when /3 = — | is the leading order term of ( l35i) non-constant in the limit t — )■ 00, and in 
this case, after also substituting back r using ( 133|) . 

J^{T;t;aL,aR) ^ e +0 - . (36) 

Finally summing the non- neglected terms yields ( 129|) for a^^aR > 0. 

When aL,ciR > 0, ( l29i) is the PDF of a Gaussian with mean ^^r^^j ^ind standard deviation 

aj^au ' -^^S- E III this case the direction of the drift is towards the origin for both positive and 
negative x. The fraction of time spent in [0, 00), approaches a^+aR ^-s t — ?■ 00. This is consistent 
with the fact that the area under the steady-state density of x(t) over [0, 00) is ^^^^^ , [2]. 



4 An application to stochastically perturbed sliding mo- 
tion 

Filippov systems are discontinuous vector fields used to model physical systems involving a dis- 
continuity or abrupt change, such as mechanical systems with impacts and relay control systems 
[T8l [19] . Stable sliding motion refers to dynamics that are constrained to a discontinuity surface 
due to a vector field that points towards the surface from both sides. 

Consider the following two-dimensional stochastic differential equation. 



dxi 
dx2 



hL 

-ClR 
bR 



xi <0 



xi > 



dt + y/e 



dWi{t) 

KdW2it) 



(37) 



where Wi{t) and W2(t) are independent standard Brownian motions, a^.aR > 0, bL,bR G M and 
e,K > 0. When e = 0, ( l37l) is Filippov and trajectories slide along xi = because ai^cLR > 0. 
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For small e > 0, f l37|) describes stochastically perturbed sliding motion and X2{t) represents 
displacement along the discontinuity surface. We can write 

X2{t) = X2{0) + bLt + {Br - Bl) [ X[0,oo) (^l (s) ) rfs + ^2 W ■ (38) 

Jo 

Consequently the PDF of X2{t) is a convolution involving the PDF of the positive occupation 
time of Xi{t) with the PDF of W2{t). Specifically, by using Theorem [1] and appropriately accom- 
modating for e, we have 



oo 



„2 



n V.^ I - I , X2 - X2i0) - bit - W t Xi(0) \ 1 

Prob t G S = / — — / pi — ; -; , a^, — == e^-* 

(39) 

By (1551) . and the asymptotic calculations of ^ as t — )■ oo, 

E(x2(t)) ^ X2(0) + 6z.t + (6«-6L)^^=X2(0) + ^^^^^^t^t, (40) 

aL + aR 0,1 + clr 

Var(x2(t)) ^ fy^^^^ + '^V^' (41) 
which agrees with results of [2J that were obtained from the PDF of ([I]). 



5 Discussion 

This paper gives for the first time the PDF of the positive occupation time of Brownian motion 
with two- valued drift. The PDF, p, is given by the complicated expression (E]). As t — )■ oo, p 
approaches the simpler form fl2^ . In the case that the piecewise drift points towards the origin 
on both sides, (129|) is Gaussian. The results help us understand the system (137|) that describes 
stochastically perturbed sliding motion of a piecewise-smooth vector field. An outstanding prob- 
lem is the computation of the two-dimensional joint PDF of x{t) and r which would enable us 
to obtain the two-dimensional transitional PDF of (I57|) . 
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